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1  Summary 


Accurate  estimation  of  the  bistatic  radar  cross  section  of  large  vehicles  is  important  for  the  design 
of  advanced  radar  systems  that  employ  spatially-distributed  transmitters  and  receivers.  These 
advanced  systems  can  enhance  radar  performance  through  the  use  of  multiple  input/multiple 
output  (MIMO)  operation  while  reducing  the  vulnerability  to  countermeasures  such  as  jamming 
and  anti-radiation  weapons.  However,  the  measurement  of  bistatic  radar  cross  sections  of  large 
vehicles  with  far-held  techniques  can  be  time  consuming,  expensive,  and  requires  physically  large 
and  less  secure  measurement  facilities.  The  measurements  must  be  made  over  a  large  set  of 
angles,  frequencies,  and  polarizations  for  each  orientation  of  the  vehicle  with  respect  to  potential 
transmit  sources.  Alternatively,  near-held  scattering  measurements  permit  estimation  of  bistatic 
scattering  using  smaller  and  more  secure  facilities.  The  use  of  advanced,  electromagnetic 
processing  enables  the  transformation  of  the  multiple,  near-held  measurements  to  the  required 
complete  set  of  far-held  cross  sections.  Still,  near-held  measurements  are  time  consuming  and 
therefore  expensive.  The  objective  of  this  work  is  to  investigate  the  feasibility  of  signihcantly 
reducing  the  time  to  make  near-held  scattering  measurements  through  the  use  of  an  array  of 
near-held  probes  instead  of  the  conventional  single  probe. 


The  use  of  an  array  of  near-held  probes  introduces  technical  problems.  Coupling  between  the 
probe  array  and  the  scattering  body  may  introduce  errors  in  the  near-held  measurements. 
Currents  induced  on  the  probe  array  and  its  support  structure  will  reradiate  and  change  the  plane 
wave-induced  currents  on  the  scattering  body.  Further,  mutual  coupling  between  the  probe  array 
elements  in  the  presence  of  the  scattering  body  may  change  the  near  held  measured  at  each  probe 
array  element  from  that  measured  by  an  isolated  probe.  In  this  work,  we  analyze  these  technical 
problems  and  describe  an  approach  to  compensate  for  the  errors  introduced  by  the  use  of  the 
probe  array. 


Based  on  an  impedance  model  for  the  interactions  between  the  probe  array  and  the  scattering 
body,  we  have  developed  expressions  for  the  measured  voltages  at  each  of  the  transmission  line 
feeds  to  the  probe  array  elements  in  terms  of  (1)  the  open-circuit  voltages  induced  by  the  plane 
wave-excited  body  at  each  probe  array  element  and  (2)  the  in  situ  mutual  impedances  of  the 
probe  array.  The  expressions  show  a  linear  set  of  relations  to  the  open-circuit  voltages  as  would 
be  expected  by  the  superposition  principle  applied  to  this  linear  system.  The  linear  set  of 
equations  can  be  inverted  to  yield  the  open-circuit,  plane  wave-induced  voltages  at  each  probe 
array  element.  We  hypothesize  that  these  voltages,  when  corrected  for  the  direct  plane  wave 
excitation,  are  linearly  related  to  the  required  near-field  measurements.  We  also  describe  a 
separate  set  of  vector  network  analyzer  measurements  to  yield  the  in  situ  mutual  impedances  (the 
mutual  impedances  in  the  presence  of  the  scattering  body)  that  are  required  to  invert  the 
measurement  equations. 


Simple  examples  of  the  importance  of  the  open-circuit  conditions  in  the  measurement  of  near 
fields  with  a  dipole  probe  array  are  given.  The  results  shown  in  the  examples  are  based  on  a 
method  of  moments  solution  to  the  electromagnetic  interaction  problem  using  NEC4.  In  these 
examples  the  scattering  target  is  considered  to  be  an  array  of  half-wavelength  dipoles  both  near 
and  far  from  the  probe  array.  The  target  dipoles  are  a  convenient  surrogate  for  a  more  complex 
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target  since  perturbations  in  the  induced  currents  due  to  the  presence  of  the  probe  array  are 
easily  found.  The  open-circuit  conditions  for  the  elements  of  the  probe  array  minimize  the 
dominant  currents  on  the  elements  leaving  smaller,  residual  currents  on  the  remaining 
quarter- wavelength  sections  of  the  dipoles.  Minimization  of  the  dominant  currents  in  the  probe 
array  elements  reduces  the  errors  in  the  target  dipole  currents  by  more  than  20  dB  when 
compared  to  the  free  space-induced  currents. 


The  preliminary  results  shown  here  are  encouraging  and  suggest  further  analysis  of  open  issues. 
These  issues  include  (1)  calibration  of  the  open-circuit,  plane  wave-induced  voltages  and  their 
relationship  to  the  required  near-field  measurements,  (2)  necessity  and  accuracy  requirements  for 
in  situ  near-field  impedance  measurements  to  characterize  the  probe  array  as  its  positions  changes 
with  respect  to  the  scattering  body,  and  (3)  experimental  and/or  numerical  confirmation  of  the 
theory  described  in  this  report. 


2  Introduction 

As  part  of  a  program  to  advance  the  art  of  bistatic  radar  cross  section  measurements,  we  consider 
the  use  of  an  array  of  field  probes  to  facilitate  simultaneous,  near-field  measurements  of  the 
electromagnetic  scattering  from  which  bistatic  cross  sections  are  estimated.  A  significant  number 
of  near-field  measurements  are  required  for  bistatic  cross  section  determination  and  the  use  of  an 
array  of  measurement  probes  will  substantially  reduce  the  time  for  bistatic  cross  section 
characterization.  However,  the  electromagnetic  interaction  between  an  array  of  measurement 
probes  and  the  scattering  body  may  be  significant  enough  to  modify  the  plane  wave-induced 
current  distribution  on  the  scattering  body  to  be  measured,  thus  causing  errors  in  the  estimated 
bistatic  cross  sections.  The  objective  of  this  note  is  to  investigate  an  approach  to  minimize  the 
interaction  between  the  array  probe  elements  and  the  nearby  scattering  body. 


Specifically,  we  consider  an  array  of  dipole  field  probes  located  near  a  scattering  body  when  the 
body  is  illuminated  by  a  plane  wave  source.  Both  the  plane  wave  source  and  the  body  induce 
currents  in  the  dipole  field  probes  which,  in  turn,  reradiate  to  the  scattering  body.  Since  the  field 
probes  are  located  near  the  scattering  body,  we  expect  that  the  probe  reradiated  fields  will 
disturb  the  currents  on  the  scattering  body  as  induced  by  the  plane  wave  excitation,  causing 
errors  in  the  scattered  field  measurements.  Further,  the  field  samples  at  each  probe  are  influenced 
by  the  mutual  coupling  to  other  elements  in  the  probe  array,  thus  introducing  errors  in  the  field 
samples.  The  mutual  coupling  between  elements  of  the  probe  array  depends  on  the  orientation  of 
the  array  with  respect  to  the  scattering  body  and  are  likely  to  change  as  the  array  orientation 
with  respect  to  the  scattering  body  changes. 

In  this  work  we  investigate  the  following  basic  questions: 

•  Is  it  possible  to  reduce  the  errors  in  probe  field  samples  caused  by  mutual  coupling  between 
the  array  elements? 

•  Is  it  possible  to  minimize  the  scattering  from  the  probe  array  in  the  near- field  volume  occupied 
by  the  scattering  body? 
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In  response  to  the  second  question,  we  analyze  an  approach  to  minimizing  the  dominant  mode 
currents  in  the  dipole  elements  of  the  probe  array,  thus  minimizing  the  reradiation  to  the 
scattering  body.  In  the  process  of  minimizing  the  dominant  mode  currents,  the  effects  of  the  in 
situ  mutual  coupling  between  array  elements  must  be  considered,  in  effect  resolving  the  issue 
raised  in  the  first  question. 


3  Model  for  the  Probe  Array  Located  in  the  Scatterer  Environment 

We  begin  by  viewing  the  near-field  measurement  system  as  a  linear  system  characterized  by 
mutual  impedances  between  elements  of  the  system  [1]  as  illustrated  in  Figure  1.  Here  the  plane 
wave  source  excites  the  scattering  body  and  the  measurement  background  which  is  characterized 
by  an  unknown  collection  of  virtual  scattering  sources.  The  plane  wave  source  also  excites  the 
measurement  array  consisting  of  a  collection  of  field  probes  terminated  in  sources  and/or  load 
impedances.  The  mutual  coupling  between  array  elements  is  characterized  by  the  in  situ  mutual 
impedance  array,  Za,  giving  the  voltage  induced  in  each  probe  array  element  by  the  currents  in 
the  other  array  elements  when  the  array  is  in  the  presence  of  the  scattering  body.  Note  that  we 
have  made  the  tacit  assumption  that  there  is  no  coupling  between  the  scattering  body  and  the 
plane  wave  source  and  between  the  array  and  the  plane  wave  source.  This  is  because  the 
scattering  from  the  body  and  array  is  weak  at  the  extended  range  of  the  plane  wave  source. 

The  details  of  the  scattering  source  are  unknown  and  change  with  the  orientation  of  the 
scattering  source  with  respect  to  the  measurement  array  and  the  plane  wave  source.  For  that 
reason,  we  focus  on  characterizing  the  coupling  between  the  plane  wave  source,  Zpa ,  and  the 
coupling  between  array  elements  in  the  presence  of  the  scattering  source,  Za- 


Virtual  Scattering  Sources 


Figure  1:  Impedance  Characterization  of  Probe  Array  Measurement  System 


We  characterizing  the  array-source-scatterer  in  terms  of  an  N  +  2-port  linear  network.  Here  there 
will  be  N  array  elements.  One  of  the  remaining  ports  is  the  source  of  plane  wave  excitation  of  the 
target  and  probe  array  and  the  other  “port”  is  associated  with  a  vector  of  ports  associated  with 
the  natural  mode  currents  on  the  scattering  body. 
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For  simplicity  throughout  this  paper  we  will  consider  an  N  =  3  element  array  since  the  results  for 
larger  arrays  will  be  an  obvious  extension.  Then  the  impedance  matrix  characterizing  the  three 
element  probe  array  measurement  system  is 


Ep 

— J- 

Zp,p  Zp,s  ^P,  1  2  Zp,  3 

Ip 

Es 

Z s,p  Z Sfs  Za9l  Zs2  Zs% 

Is 

Ei 

= 

Z\,p  Zps  Z\  \  Z\;2  Z 13 

II 

E'2 

Z2,p  Z2s  Z2.\  Z2,2  Z2;3 

h 

e3 

Z3,p  Z3s  Z3,1  Z3,  2  Z;i3 

h 

Here  Ep.  Es  and  Ei,  i  =  1,2,3  are  the  plane  wave  source,  scattering  body  natural  mode,  and 
dipole  voltages,  respectively  and  Ip,  Is ,  and  i  =  1,2,3  are  the  corresponding  terminal  currents. 
The  impedances,  Z.;j,  i .  j  =  1,2,3,  are  the  mutual  and  self  impedances  of  the  array  elements  with 
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the  natural  modes  of  the  scattering  body  and  the  column  matrix  ZStP  characterizes  the  coupling 
between  the  natural  modes  and  the  plane  wave  source.  Finally,  the  impedance  column  matrices 
Zs< i,  Zs> 2,  and  Zs )3  characterized  the  coupling  between  the  natural  modes  of  the  scattering  body 
and  the  elements  of  the  probe  array. 


We  focus  on  the  last  three  rows  of  this  matrix  characterizing  the  probe  dipole  voltages  and 
currents  and  the  open-circuit  voltages  induced  in  the  dipole  elements  by  the  plane  wave.  That  is, 
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Now  the  currents  in  the  scatterer  natural  modes  can  be  found  from  the  second  row  equation  of  (1) 
by  assuming  that  the  scatterer  terminals  are  loaded  with  equivalent  impedances  Zs  giving 
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Substituting  into  Equation  2  gives 
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Here  Ef ,  i  =  1,2,3  represent  the  open-circuit  voltages  induced  in  each  probe  array  element  by 
the  incident  plane  wave  source  in  the  presence  of  the  scattering  body  and  Zij,  i,j  =  1, 2,  3  are  the 
in-situ  mutual  impedances  of  the  probe  array.  In  neglecting  consideration  of  the  first  equation  we 
have  assumed  that  any  radiation  from  the  dipole  elements  in  the  presence  of  the  scattering  body 
and  the  scattering  from  the  body  itself  does  not  change  the  primary  source  of  plane  wave 
excitation,  as  noted  earlier. 

It  is  important  to  note  that  each  of  the  the  plane  wave-induced,  open-circuit  voltages,  Ef,  is  the 
superposition  of  (1)  a  voltage  due  to  the  direct  excitation  of  the  probe  array  element  by  the  plane 
wave,  ZitPIp,  and  (2)  a  voltage  due  to  the  scattering  from  the  body  when  excited  by  the  plane  wave, 
zfs(ZS)S  +  Z s)~l Z S'plp  .  In  this  report,  we  will  develop  an  approach  to  measuring  the  open-circuit 
voltages,  Ef ,  which,  when  corrected  by  removing  the  direct  excitation  term,  are  proportional  to 
the  near-zone  scattered  fields  at  the  probe  array  elements. 

Equation  3  suggests  the  N-port  representation  of  the  probe  array  measurement  system  as 
illustrated  in  Figure  2.  Here  the  ports  are  the  probe  array  element  terminals  with  voltages  Ejt  and 
currents  Ii,i  =  1,2,...,  N.  Each  port  has  the  open-circuit  voltage,  Ef ,  induced  by  the  plane 
wave  source  in  the  presence  of  the  scattering  body. 


Figure  2:  N-Port  Equivalent  Circuit  of  Measurement  System 


4  Minimizing  Probe  Array  Interactions 


We  develop  an  approach  to  actively  load  the  terminals  of  the  probe  array  to  minimize  array 
reradiation  which  can  modify  the  plane  wave  induced  currents  on  the  scattering  body.  In  the 
following  sections  we  develop  a  theoretical  basis  for  minimizing  the  currents  in  the  probe  array 
elements  and  present  an  elementary  example  to  demonstrate  performance. 


4.1  Basic  Theory 


The  objective  in  this  work  is  to  actively  load  the  terminals  of  the  probe  array  to  minimize  array 
reradiation  which  can  modify  the  plane  wave  currents  on  the  scattering  body.  To  that  end  we 
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consider  the  model  for  the  N-port  measurement  system  with  Thevenin  equivalent  source  voltages 
Esi  and  impedances  ZSi,  i  =  1,2,...,  N  at  the  probe  array  terminals  as  shown  in  Figure  3. 


— 0^- 

Passive 

N  port  Network 

• 

• 

-•“SD-p) 

Figure  3:  N-Port  Measurement  System  with  Thevenin  Equivalent  Sources  at  each  Terminal 


In  this  case,  the  simplified  network  equations  given  by  Equation  (3)  become 


Ex 

Esi  —  Zsih 

E[ 

’  ?i,i 

Z\,2 

Z\;3 

h 

e2 

= 

Es  2  —  Zs2h 

= 

E‘2 

+ 

Z2,i 

Z-2,2 

Z2,3 

h 

e3 

Es3  —  Zs3I3 

Z3}\ 

Z3,2 

z3,3 

h 

which  is  equivalent  to 


'  Es  1  -  Ef  ' 

^1,1  +  Zs  1  Z\2  Z\3 

h 

Es 2  -  E? 

= 

Z2,l  Z2,2  +  ZS2  Z2j3 

h 

i 

1 

CO 

r-  ^ 

KJ 

_ 1 

z3,l  Z3‘2  Z33  +  ZS3 

h 

(4) 


It  is  clear  from  these  equations  that  forcing  all  elements  of  the  column  matrix  on  the  left  side  to 
be  equal  to  zero  assures  that  the  currents  induced  in  each  probe  element  port  will  be  zero, 
providing  of  course  that  the  augmented  impedance  matrix  is  not  singular.  Requiring  that  the 
column  matrix  to  the  left  be  zero  is  equivalent  to  requiring  that  the  open-circuit  Thevenin 
equivalent  source  voltage  at  each  port  be  equal  to  the  open-circuit  voltage  in  each  probe  array 
element  induced  by  the  plane  wave  excitation  in  the  presence  of  the  scattering  body. 


4.2  Elementary  Example 


To  illustrate  this  point  we  consider  the  following  elementary  numerical  example.  Here  a  single 
measurement  dipole  is  located  near  a  collection  of  passive,  target  dipoles  which  serve  to  represent 
the  measurement  body.  The  measurement  dipole  and  target  dipoles  are  illuminated  by  a  plane 
wave.  The  use  of  dipoles  is  both  convenient  and  practical.  The  electromagnetic  environment  will 
be  characterized  by  the  Numerical  Electromagnetic  Code  (NEC4)  which  conveniently  determines 
the  interaction  between  collections  of  dipoles  which  may  be  loaded  and  excited  by  local  sources  as 
well  as  a  plane  wave  source.  Further,  the  dipole  is  a  practical  measurement  probe  used  in  some 
antenna  and  cross  section  ranges.  Also,  dipoles  are  a  convenient  surrogate  for  the  scattering  body 
since,  with  the  NEC  code,  perturbations  of  the  currents  on  the  scattering  dipoles  can  be 
determined  as  a  measure  of  the  magnitude  of  the  probe-scatterer  interactions.  In  this  example, 
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the  measurement  dipole  will  be  fed  by  a  source  at  its  center.  We  will  determine  the  required 
source  voltage  to  cancel  the  induced  current  on  the  measurement  dipole  and  then  examine  the 
current  perturbations  on  the  scattering  dipoles. 


Specifically,  the  single  measurement  probe  consists  of  a  thin  dipole  of  length  0.48A,  where  A 
denotes  the  free-space  wavelength  of  the  measurement  plane  wave.  The  probe  dipole  can  be 
excited  by  a  voltage  source  at  its  center.  The  voltage  source  is  assumed  to  be  coherent  with 
respect  to  the  plane-wave  illumination.  That  is,  the  frequency  of  the  voltage  source  is  that  of  the 
plane  wave  and  its  phase  is  constant  during  measurements.  As  illustrated  in  Figure  4,  two  target 
dipoles  are  located  adjacent  to  the  measurement  probe.  Both  are  half-wavelength  dipoles.  One  is 
located  5A  from  the  probe  and  the  other  is  located  at  35A  away  from  the  probe  and  offset  by  5A. 
The  collection  of  probe-target  dipoles  is  excited  by  a  plane  wave  with  an  E  field  strength  of  lv/rn. 


Figure  4:  Geometry  for  Simple  Example  of  Single  Probe  Element  and  Target  Dipoles 

We  use  NEC4  to  determine  relevant  properties  of  the  probe  and  target  dipole  current 
distributions.  First,  we  determine  the  short-circuit  current  due  to  plane  wave  excitation  at  the 
center  of  the  probe  dipole  when  the  target  dipoles  are  present.  This  is  done  by  requiring  that  the 
probe  dipole  voltage  source  be  zero  volts.  Second,  we  determine  the  current  excited  by  the 
voltage  source  at  the  center  of  the  probe  dipole  when  the  target  dipoles  are  present  and  there  is 
no  plane-wave  excitation.  Finally,  using  superposition  in  this  linear  system  with  both  plane  wave 
and  voltage  source  excitations,  the  voltage  is  adjusted  to  create  a  current  that  cancels  that 
induced  by  the  plane  wave.  The  effectiveness  of  the  cancelation  is  determined  by  comparing  the 
currents  in  the  target  dipoles  with  those  induced  in  the  dipoles  when  the  probe  is  removed. 
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Meaurement  Dipole  Currents 
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- with  Voltage  Source 

-  -  -  with  Cancellation  Source 


Figure  5:  Current  Distributions  on  Measurement  Dipole  in  Simple  Example 


Figure  6:  Current  Distributions  in  Target  Dipoles  and  Current  Errors  Due  to  Presence  of  Measure¬ 
ment  Dipole  with  Canceled  Current 

Figures  5  and  6  show  the  results.  Figure  5  shows  the  measurement  dipole  currents  with  plane 
wave,  voltage,  and  cancelation  sources.  The  voltage  source  is  set  at  1  volt.  Note  that  the  current 
is  effectively  canceled  at  the  center  of  the  measurement  dipole  where  the  cancelation  voltage  is 
applied.  However,  a  residual  current  remains  which  is  approximately  10  dB  below  the  current 
induced  by  the  incident  plane  wave.  While  this  does  not  appear  to  be  a  substantial  cancelation  of 
the  measurement  dipole  current,  the  real  measure  of  effectiveness  is  the  perturbation  in  currents 
induced  in  the  target  dipoles.  These  results  are  shown  in  Figure  6. 


The  target  dipole  currents  with  no  measurement  dipole  present  and  with  the  canceled  current  on 
the  measurement  dipole  are  shown  in  the  upper  row.  Since  the  difference  between  the  currents  is 


too  small  to  note  in  these  graphs,  plots  of  the  magnitude  of  the  error  in  the  complex  currents  are 
presented  in  the  second  row.  For  the  near  dipole  (5A  away  from  the  measurement  dipole),  the 
error  is  less  than  -50  dB  below  the  magnitude  of  the  current  on  the  dipole.  For  the  distant  dipole 
(35 A  away  from  the  measurement  dipole),  the  error  is  less  than  -65  dB  below  the  target  dipole 
current  magnitude. 


We  conclude  that  even  a  modest  reduction  in  the  measurement  dipole  currents  results  in  very 
small  errors  in  the  target  dipole  currents  when  compared  to  currents  induced  in  the  free-space 
environment  alone.  This  is  an  encouraging  illustration  of  the  approach  described  here  and  forms 
the  basis  for  the  following  analysis  of  a  more  realistic  measurement  system. 


4.3  Open-Circuit  Induced  Voltages  on  Probe  Array  Elements 


It  is  clear  from  the  previous  discussion  of  Equation  4  that  the  plane  wave-induced  open-circuit 
voltage  at  each  terminal  of  the  probe  array  must  be  measured.  In  the  following  we  suggest  an 
approach  to  making  those  measurements. 


We  begin  by  adopting  a  more  realistic  model  for  the  N-port  network  representing  the 
measurement  array  (Figure  3)  where  now  each  network  terminal  is  connected  to  its  Thevenin 
equivalent  source  through  a  transmission  line.  This  model  is  shown  in  Figure  7.  To  understand 
the  process  by  which  we  measure  the  plane  wave-induced  open-circuit  voltages,  we  develop  the 
Thevenin  equivalent  circuit  for  the  N-port  probe  array  network  at  each  network  port  with  the 
transmission  lines  and  loads  at  the  other  ports.  This  is  illustrated  in  Figure  8  where  the  arrow 
indicates  the  point  at  which  the  Thevenin  equivalent  circuit  is  to  be  developed  at  network 
terminal  1  and  the  other  transmission  lines  are  shown  loaded  by  their  Thevenin  equivalent  source 
impedances. 


Figure  7:  N-Port  Equivalent  Circuit  of  Measurement  System  with  Transmission  Line  and  Sources 
at  Each  Port 
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Figure  8:  N-Port  Equivalent  Circuit  for  Finding  Thevenin  Equivalent  Circuit 


We  develop  the  network  equivalent  circuit  at  port  1  assuming  a  three-port  network.  From 
Equation  (3)  we  can  write  the  following  equations  which  determine  the  equivalent  open-circuit 
voltage,  E°c,  of  the  Thevenin  equivalent  circuit 


Ei 

E\c 

E[ 

Z\;z 

z  1 .3 

0 

e2 

= 

—Z2I2 

= 

E? 

+ 

Z2,l 

•^2,2 

■^2,3 

h 

E3 

-Z3I3 

E3 

■Z’a,] 

•^3,2 

-^3,3 

h 

where  by  setting  I±  =  0  we  find  the  open-circuit  voltage  at  port  1.  Note  that  the  other  network 
ports  are  loaded  with  impedances  Z2  and  Z3  and  the  negative  signs  account  for  the  convention 
that  the  port  currents  are  assumed  to  be  positive  when  they  enter  the  port  terminal.  The  load 
impedances  are  easily  found  from  the  transmission  line  equations  [1]  at  each  port  as 


Zi  —  Z a 


1  +  rie-j2kLi 

1  -  Ti.e-i2kLi  ’ 


(5) 


where 


r,  = 


Zsi  ~  Z0 
Zsi  +  Z0 


and  Z0  =  the  characteristic  impedance  of  the  transmission  line,  Lj 
transmission  line  at  port  i  and  k  = 


the  length  of  the 


It  follows  from  the  first  of  the  set  of  equations  that 


E 


OC 

1 


E?  +  Z\p.  Z 1  3 


Ei  —  Zip  Zi.3 


h 

h 


J 

•^2,2  +  Z2 

^2,3 

-1 

'  E^  ' 

Z3,2 

Z3,3  +  Z3 

.  Ez  . 

(6) 


Here  the  second  and  third  of  the  set  of  equations  have  been  solved  for  the  currents  I2  and  I3  and 
substituted  into  the  first  equation. 
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To  develop  the  Thevenin  equivalent  impedance,  we  determine  the  short-circuit  current  at  port  1, 
I[c,  since  the  Thevenin  equivalent  impedance  is  the  ratio  Zyth  =  E°c/I(c.  Again,  referring  to 
Equation  3  we  can  write 


Ei 

0 

E[ 

'  2U 

2i,2 

2i,3 

JSC 

— il 

e2 

— 

—Z2I2 

= 

E'2 

+ 

22,1 

22,2 

22,3 

h 

E3 

—Z3I3 

Ei 

23,1 

23,2 

23,3 

h 

The  short-circuit  current  at  port  1  is  determined  from  the  first  of  the  set  of  equations  by  requiring 
that  Ei  =  0.  Note  that  the  current  I\  =  —Ifc  where  the  minus  sign  is  included  to  maintain  the 
current  convention  for  the  Thevenin  equivalent  circuit.  From  the  second  and  third  of  these 
equations  we  find 


Z2.1 

Z3.1 


JSC  _ 

~ 


Ei 

Ei 


+ 


Zi2  _  -Z-2,3 

•Z3.2  ^3,3  +  Zj, 


h 
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Then  from  the  first  equation  we  find 
E\  =  2i,iT(c  —  [  Zl,2  2^1,3 


=  I  Z\  \  - 


Zi,ilfc  —  |  Zi,2  Z1>3 
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Using  the  expression  for  E°c  in  Equation  6,  this  result  becomes 


ur  =  1  Z\  \  - 


2i,2  2i,3 


22,2^+  Z2  ~  22,3 
23,2  23  3  +  23 


-1 


22, ]_ 

23,1 
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where  the  Thevenin  equivalent  impedance  at  port  1  is 

17  _  TPoc  /  Tsc  _  V 

^lTh  ~  Hi  ~  ^1,1  — 


ry  ry 

22,2  +  Z2  22,3 

-1 

22,1 

^1,2  ^1,3 

23,2  23,3  +  23 

23,1 

(7) 


Note  that  the  Thevenin  equivalent  circuit  for  the  loaded  N-port  (Equations  (5,  6,  and  7)  require 
that  we  know  the  in  situ  self  and  mutual  impedances  of  the  probe  array  as  well  as  the  plane  wave 
induced  open-circuit  voltages  at  each  of  the  probe  array  ports. 


To  complete  the  analysis  we  transform  the  Thevenin  equivalent  source  through  the  transmission 
line  at  port  1  to  the  source  end  of  the  transmission  line.  Here  we  require  that  the  source  end  of 
the  transmission  line  be  loaded  with  the  source  impedance  Zs\  and  we  set  the  source  voltage 
Es  1  =  0  volts  as  illustrated  in  Figure  9.  With  standard  transmission  line  analysis  [2]  we  find  that 
the  measured  voltage  at  the  source  end  of  the  transmission  line  is 


E\m  — 


EfZoil  +  T,) 


(ZiTh  +  20)(1  -  YThVie-i2kL^) 


,-jkL-i 
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Figure  9:  N-Port  Equivalent  Circuit  with  Measurement  Load 


where 


r™  — 


Z\Th  —  Z0 


ZlTh  +  Z0 

Substituting  for  E°c  from  Equation  6  gives  the  measured  voltage  at  source  port  1  as 
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In  this  equation,  A  denotes  the  determinant  of  the  reduced  impedance  matrix: 


(8) 


A  = 


Z2,2~+  Z2  ~  Z2,3 
Z3,2  Z33  +  Z3 


It  is  important  to  note  in  Equation  8  that  the  measured  voltage  at  source  port  1  is  a  linear 
combination  of  the  open-circuit,  plane  wave  induced  voltages  at  the  network  ports.  Of  course  this 
is  to  be  expected  given  that  the  measurement  system  is  linear  and  the  measured  response  must 
then  be  a  linear  combination  of  the  system  excitation  voltages.  By  extension,  the  measured 
voltages  at  the  other  measurement  ports,  E2m  and  E3m,  are  linear  combinations  of  the  same 
open-circuit,  plane  wave  induced  voltages.  Thus,  in  the  case  of  the  three  port  system,  we  have 
three  linear  equations  which  can  be  inverted  to  give  the  three  required  plane  wave  induced 
voltages.  Note  however,  that  this  inversion  requires  that  we  measure  independently  the  in  situ 
impedance  matrix  given  in  Equation  3  when  the  plane  wave  excitation  is  turned  off.  We  will 
describe  an  approach  to  measuring  these  impedances  later  in  this  paper. 


Before  we  leave  this  development,  we  determine  the  source  voltage  required  to  cancel  the  current 
at  the  probe  network  terminals  when  the  probe  array  is  excited  by  a  plane  wave.  Consider  Figure 
10  where  measurement  port  1  is  driven  by  an  active  source  with  equivalent  Thevenin  open-circuit 
voltage  Es\  and  impedance  Zs  \ .  The  probe  array  plane  wave  excitation  is  characterized  by  the 
induced  open-circuit  voltages  Ef,  i  =  1,  2, . . .  i  N. 
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If  we  transform  the  Thevenin  equivalent  source  to  the  probe  array  network  port  1,  a  new 
equivalent  source  is  formed  at  the  probe  array  port  with  the  open-circuit  voltage  and  equivalent 
impedance  [3] 

=  2  EslZ0e~^ 

lseq  (zsl  +  z0)(i-Tlseqy 

and 

_  Z\s  +  j Z0tan(kLi)  _  1  +  T  lseq 

lseq  °  ZQ  + jZistan(kLi)  °l-YUeq' 

where  =  Tlse~j2kLl  and  Tls  =  . 


Figure  10:  N-Port  Equivalent  Circuit  with  Active  Load 


Figure  11  illustrates  the  equivalent  circuit  at  probe  array  port  1  where  the  equivalent  circuit  of 
the  loaded  probe  array  with  plane  wave  illumination  is  shown  to  the  left  and  the  equivalent 
circuit  for  the  source  at  measurement  port  1  is  shown  to  the  right.  It  is  clear  from  this  illustration 
that  to  drive  the  port  current  to  zero  requires  that  E°c  =  E\seq  or 
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giving  the  required  voltage  a  measurement  port  1  as 
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Figure  11:  Load  Voltage  to  Cancel  Probe  Current 
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Again  we  note  that  the  cancelation  voltage  Es\  requires  knowledge  of  not  only  the  in  situ  self  and 
mutual  impedances  of  the  probe  array  but  also  the  plane  wave  induced  open-circuit  voltages  at 
each  of  the  probe  array  ports. 


Recall  that  the  objective  of  the  probe  array  is  the  simultaneous  measurement  of  the 
near-scattered  field  at  each  probe  array  element.  It  seems  that  the  solution  for  the  plane  wave 
induced  open-circuit  voltages,  Ef ,  i  =  1,2, ... ,  N ,  using  Equation  8  and  the  measured  voltages 
Eim ,  i  =  1,2, ,  N  would  satisfy  this  objective,  once  the  open-circuit  voltages  have  been 
corrected  for  the  term  characterizing  the  direct  plane  wave  excitation.  Note  that  in  this  process 
we  have  determined  the  open-circuit  voltages  mathematically  without  physically  inducing  the  zero 
current  at  each  probe  network  port  as  implied  in  Equation  9. 


4.4  Determination  of  Probe  Array  In  Situ  Impedances 


Finally,  we  develop  an  approach  to  measure  the  in  situ  self  and  mutual  impedances  of  the  probe 
array  which,  as  noted  earlier,  are  required  to  determine  the  plane  wave  induced  open-circuit  probe 
array  voltages.  We  adopt  the  basic  measurement  technique  used  in  modern  vector  network 
analyzers.  Consider  the  probe  array  network  shown  in  Figure  12  with  equivalent  sources  and 
transmission  line  feeds.  We  assume  that  there  is  no  plane  wave  excitation  for  these  measurements 
and  that  the  source  impedance  matches  the  transmission  line  characteristic  impedance,  ZQ.  The 
signals  on  each  transmission  line  are  characterized  by  their  incident  and  reflected  voltages,  a™ 
and  6™,  i  =  1,2, ,  N  respectively. 
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Figure  12:  N-Port  Equivalent  Circuit  for  Finding  Array  Mutual  and  Self  Impedances 

The  linear  relationship  between  the  incident  and  reflected  voltages  is  characterized  by  the 
scattering  matrix  for  the  network  and  transmission  lines  as 
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Here  the  superscript  m  denotes  that  these  parameters  refer  to  the  measurement  ports  of  the 
probe  array  network  and  transmission  line  system.  Each  of  the  scattering  parameters  can  be 
determined  from  direct  measurement  of  the  incident  and  scattered  voltages  as,  for  example, 

hm 

m  ^2  | 

S23  =  |o™=a™=0  • 

a3 

The  constraint  a™  =  a™  =  0  implies  that  the  sources  at  measurement  ports  1  and  2  are  removed 
and  all  ports  are  terminated  with  the  characteristic  impedance  of  their  transmission  lines. 


We  require  the  scattering  matrix  and  associated  mutual  impedance  matrix  for  the  probe  array.  To 
that  end,  we  de-embed  the  incident  and  reflected  voltages  at  the  measurement  ports  to  the  probe 
array  ports  yielding  af  and  bf,i  =  l,2,...,N.  That  is, 

of  =  e-jkLia ™  (11) 


and 


bf  =  ejkL'b™, 


(12) 


for  all  *  =  1,2,...,  N.  Of  course  the  incident  and  reflected  voltages  at  the  probe  array  ports  yield 
the  scattering  matrix  for  the  probe  array  as 
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If  we  use  Equations  11  and  12  with  Equation  10  we  can  write 
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providing  an  expression  for  the  scattering  matrix  at  the  probe  array  ports  in  terms  of  the 
measured  scattering  matrix  as 
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Finally,  we  can  convert  the  scattering  matrix  at  the  probe  array  ports  to  the  array  mutual 
impedance  matrix.  Recognizing  that  at  each  probe  array  port  the  transmission  line  voltages  and 
currents  are  related  to  the  incident  and  scattered  wave  complex  amplitudes  by 

and 

n  =  -  (>?)• 

In  matrix  form 

-  a 

Ea  =  / Zo(aa  +  ba)  =  ~ZT  =Z=(aa-  ba). 

V 

-d 

This  matrix  equation  can  be  solved  for  the  probe  array  scattering  matrix  S  using  the  identity 
matrix  Id  as 


Z0Id{da  +  ba)  =  z\da  -  ba), 

or 

{z  +  z0Td)ba  =  (z  -  z0Jd)da, 

-CL 

so  that  the  scattering  matrix  S  is  given  in  terms  of  the  array  mutual  impedance  matrix  by 

ba  =  sV  =  (za  +  z07d)-\za  -  Z0Td)da. 

This  result  can  be  inverted  to  yield  the  probe  array  mutual  impedance  matrix  as  follows: 

Sa  =  (. Za  +  Z0Tci)-l(Za  -  Z0TcL) 

or  _  _ 

{Z  +  Z 0Id) S'  =  {Z  -  Zjd) 

giving 

2a{7d-^a)  =  z0(fd+Ea) 

so  that  _  _  _  _ 

Z  =  Z0(Id-  5°)_1(Jd  +  50).  (15) 

We  see  that  measurements  of  the  scattering  parameters  at  the  measurement  ports  can  be 
converted  to  the  scattering  parameters  at  the  probe  array  network  ports  using  Equation  14  and 
those  scattering  parameters  can  be  converted  to  the  probe  array  mutual  impedance  matrix  using 
Equation  15. 
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5  Example  Results  from  Numerical  Analysis 


We  conclude  with  an  elementary  example  consisting  of  a  three  element  probe  array.  This  example 
is  included  to  demonstrate  performance  when  the  probe  array  element  currents  induced  by  the 
plane  wave  and  body  scattering  are  canceled  by  probe  array  element  sources.  The  measure  of 
performance  is  the  error  in  the  target  dipole  currents  produced  by  the  canceled  probe  element 
sources  when  compared  to  the  free  space  currents  induced  when  the  probe  array  is  absent. 


As  illustrated  in  Figure  13,  the  probe  array  consists  of  three  dipoles  of  length  0.48A  and 
inter-element  spacing  of  A/2.  Three  surrogate  target  dipoles  of  length  A/2  are  located  5 A  away 
from  the  array  as  illustrated  and  one  target  dipole  is  35A  from  the  probe  array  and  offset  by  5A. 
The  probe  array  and  target  dipoles  are  illuminated  by  a  plane  wave,  initially  broad  side  to  the 
probe  array  and  then  75°  from  the  broadside.  Note  that  in  this  example  the  probe  array  elements 
are  driven  directly  by  voltage  sources  rather  than  via  transmission  lines. 


Figure  13:  Geometry  for  Three  Element  Probe  Array  and  Target  Dipoles 


Analysis  of  this  configuration  was  performed  using  NEC4  to  account  for  the  electromagnetic 
interaction  between  probe  array  elements  and  the  target  dipoles  as  well  as  the  in  situ  mutual 
impedances  between  the  probe  array  elements.  Initially  the  in  situ  mutual  impedances  are 
determined  using  NEC4  by  driving  each  element  separately  with  its  voltage  source  and 
calculating  the  short-circuit  current  in  the  probe  array  elements,  in  effect  determining  the  in  situ 
admittance  matrix  of  the  array.  The  array  mutual  impedance  matrix  is  the  inverse  of  the 
admittance  matrix.  For  these  impedance  calculations,  no  plane  wave  excitation  is  applied.  Next 
the  probe  array  element  voltage  sources  are  adjusted  using  the  impedance  matrix  to  cancel  the 
currents  induced  in  the  probe  array  elements  by  the  plane  wave  excitation.  In  this  state,  the 
currents  in  the  target  dipoles  are  calculated  and  compared  with  those  induced  by  the  plane  wave 
alone  without  the  probe  array  present. 
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Results  of  this  analysis  are  shown  in  Figures  14  and  15  for  broadside  plane  wave  excitation  and 
Figures  16  and  17  for  plane  wave  excitation  which  is  75°  from  broadside  illumination.  Figure  14 
shows  the  currents  in  the  short-circuited  probe  array  elements  when  excited  by  the  broadside 
plane  wave  and  the  currents  in  the  probe  array  elements  with  the  cancelation  voltage  applied. 

The  currents  on  the  end  elements  are  identical  because  of  symmetry  and  so  only  four  curves  are 
shown.  It  is  clear  that  complete  current  cancelation  occurs  only  at  the  port  where  the  voltage  is 
applied.  Residual  currents  occur  on  the  quarter  wavelength  ends  of  the  probe  element  dipoles  and 
these  residual  currents  are  approximately  20  dB  below  the  currents  on  the  short-circuit  dipoles.  It 
is  clear  once  again  that  complete  current  cancelation  with  only  a  single  load  point  is  not  possible 
leaving  residual  dipole  currents  that  can  reradiate  and  disturb  the  plane  wave  induced  currents  on 
the  target  dipoles. 


Meaurement  Dipole  Currents 


End  Element  with  Plane  Wave  Excitation 
™  ™  ™  Center  Element  with  Plane  Wave  Excitation 
End  Element  with  Plane  Wave  Excitation 
End  Element  with  Cancellation 
'  ™  '  '  Center  Element  with  Cancellation 

End  Element  with  Cancellation 


Figure  14:  Short-Circuit  and  Canceled  Currents  on  Probe  Array  Dipoles  -Broadside  Plane  Wave 
Excitation 


The  target  dipole  currents  are  shown  in  Figure  15.  The  current  distributions  on  each  of  the  target 
dipoles  shown  in  Figure  13  is  plotted  with  the  measurement  array  absent  and  with  the  canceled 
currents  on  the  measurement  array.  The  current  distribution  when  the  measurement  array  is 
absent  is  considered  the  reference  since  no  perturbation  of  the  plane  wave  induced  currents  is 
possible  in  this  case.  In  Figure  15  it  is  impossible  to  distinguish  differences  between  these  current 
distributions  and  so  we  have  plotted  the  magnitude  of  the  difference  in  the  complex  currents  in 
two  cases.  These  are  (1)  the  magnitude  of  the  error  between  the  plane  wave  induced  currents  with 
measurement  array  absent  and  with  the  short-circuited  array  present  and  (2)  the  magnitude  of 
the  error  between  the  plane  wave  induced  currents  with  the  measurement  array  absent  and  with 
the  canceled  current  array  present.  Figure  15  shows  that  canceling  the  probe  currents  reduces  the 
error  in  the  target  dipole  currents  by  20  dB  or  more.  The  magnitude  of  the  error  in  the  induced 
currents  with  the  probe  array  present  and  currents  canceled  is  between  -40  dB  and  -60  dB  with 
respect  to  the  free  space,  plane  wave  induced  currents. 

Figures  16  and  17  show  similar  results  when  the  plane  wave  excitation  is  75°  from  broadside. 

Once  again,  the  probe  element  currents  are  canceled  well  only  at  the  load  point  leaving  residual 
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currents  to  reradiate  and  modify  the  free  space  plane  wave  induced  currents  on  the  target  dipoles 
(Figure  16).  In  Figure  17  we  see  that  the  errors  introduced  by  the  probe  array  with  canceled 
currents  are  smaller  by  approximately  10  dB  in  the  center,  near-target  dipole  and  far-target  dipole 
when  compared  to  the  broadside  illumination  case.  Also  the  errors  in  the  displaced  near-target 
dipoles  (dipoles  2  and  3)  are  larger  by  20  dB  when  compared  to  the  broadside  illumination  but 
these  errors  are  still  below  40  dB  when  compared  to  the  free  space,  plane  wave  induced  currents. 


Near  Target  Dipole  1  Currents,  s  Near  Target  Dipole  2  Currents 


- No  Measurement  Array 

•  with  Cancelled  Measurement  Array  Currents 
—  —  1  Error  with  Measurement  Array  and  Cancelled  Currents 
—  —  —  Error  with  Short  Circuit  Measurement  Array 


Figure  15:  Currents  on  Target  Dipoles  and  Errors  when  Compared  to  Short-Circuit  and  Canceled 
Current  Measurement  Array  -  Broadside  Plane  Wave  Excitation 


Meaurement  Dipole  Currents 


- End  Element  with  Plane  Wave  Excitation 

—  ■  —  ■  Center  Element  with  Plane  Wave  Excitation 

- End  Element  with  Plane  Wave  Excitation 

■  End  Element  with  Cancellation 
—  —  —  Center  Element  with  Cancellation 
■'  End  Element  with  Cancellation 


Figure  16:  Short-Circuit  and  Canceled  Currents  on  Probe  Array  Dipoles  -  Plane  Wave  Excitation 
at  75°  from  Broadside 
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Near  Target  Dipole  1  Currents, 


Near  Target  Dipole  2  Currents 


- No  Measurement  Array 

•  with  Cancelled  Measurement  Array  Currents 

Error  with  Measurement  Array  and  Cancelled  Currents 
—  —  —  Error  with  Short  Circuit  Measurement  Array 


Figure  17:  Currents  on  Target  Dipoles  and  Errors  when  Compared  to  Short-  Circuit  and  Canceled 
Current  Measurement  Array  -  Plane  Wave  Excitation  at  75°  from  Broadside 


6  Conclusions  and  Recommendations 


We  have  studied  an  active  loading  approach  to  minimize  the  electromagnetic  interactions  between 
an  array  of  field  probes  and  a  nearby  scattering  body.  The  loading  minimizes  the  currents  on  the 
array  of  field  probes  while  simultaneously  measuring  the  open-circuit  voltages  induced  in  the 
probes.  We  hypothesize  that  these  open-circuit  voltages,  when  corrected  for  the  direct  plane  wave 
excitation,  are  proportional  to  the  near-zone  scattered  field  of  the  body  at  the  probe  locations. 
These  near-zone  field  measurements  are  useful  for  efficiently  synthesizing  the  far-zone  bistatic 
scattering  from  the  body. 


Numerical  studies  of  simple  examples  show  that  the  currents  on  the  approximately  half 
wavelength  probe  array  elements  are  canceled  effectively  only  at  the  single  load  point,  leaving 
residual  currents  to  reradiate  and  modify  the  plane  wave  induced  currents  on  the  scattering  body. 
Nevertheless,  examination  of  the  currents  on  surrogate  target  dipoles  indicates  that  substantial 
reductions  in  the  errors  in  the  induced  currents  results  from  the  single  point  loading  of  each 
element  in  a  three  element  probe  array. 


A  theory  for  the  active  loading  was  developed  for  a  general  case  of  the  probe  array  fed  by  active 
loads  through  transmission  lines  to  each  element  of  the  probe  array.  Measurement  of  the  plane 
wave  induced  voltages  at  each  of  the  load  ports  provides  sufficient  data  to  evaluate  the  plane  wave 
induced  open-circuit  voltages  at  each  of  the  probe  array  ports.  We  hypothesize  that  these 
open-circuit  voltages,  when  corrected  for  the  direct  plane  wave  excitation,  are  proportional  to  the 
required  near-field  samples.  However,  inversion  of  the  measured  data  requires  additional 
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measurements  of  the  in  situ  mutual  and  self  impedances  of  the  probe  array.  Here  we  use  the 
terminology  “in  situ”  to  connote  the  requirement  that  the  impedances  be  measured  when  the 
probe  array  in  located  near  to  the  scattering  body,  and  not  in  free  space. 


While  the  theory  and  simple  numerical  results  are  encouraging,  further  analysis  is  required  of 
several  issues.  These  are: 


•  Calibration  of  the  open-circuit,  plane  wave-induced  voltages.  As  indicated,  we  have 
hypothesized  that  the  open-circuit  voltages  induced  in  the  probe  array  elements  by  the  plane 
wave  and  the  scattering  body,  when  corrected  for  the  direct  plane  wave  excitation,  are 
proportional  to  the  required  near-field  samples  of  the  body  scattering.  This  hypothesis  needs  to 
be  investigated  to  insure  that  the  voltages  in  the  probe  array  elements  are  a  close  approximation 
to  the  voltages  induced  in  a  single  probe  dipole  at  the  same  location.  Further,  calibration  of  the 
probe  array  elements  (an  effective  dipole  length)  needs  to  be  developed  to  convert  the 
open-circuit  voltages  to  the  corresponding  electric  field  samples. 


•  Measurement  of  the  in  situ  impedances  of  the  probe  array.  The  self  and  mutual  impedances 
of  the  probe  array  are  a  function  of  the  geometry  of  the  probe  array  itself  as  well  as  the 
orientation  of  the  array  with  respect  to  the  nearby  scattering  body.  Since  the  probe  array  will  be 
moved  with  respect  to  the  scattering  body  to  accumulate  the  near-scattered  field  measurements, 
it  may  be  necessary  to  measure  the  self  and  mutual  impedances  at  each  probe  array  position. 
Since  these  measurements  will  extend  the  overall  measurement  time  and  the  objective  of  the 
probe  array  is  to  reduce  this  measurement  time,  we  need  to  examine  whether  and  how  frequently 
the  mutual  impedance  measurements  must  be  made.  For  example,  can  we  use  the  impedances 
measured  in  free  space  and  still  achieve  sufficient  accuracy  in  the  induced,  open-circuit  voltages? 
If  not,  how  often  must  the  impedance  measurements  be  updated  as  the  probe  array  moves  near 
the  scattering  body?  With  what  precision  do  the  impedances  need  to  be  measured? 

•  Experimental  and/or  numerical  confirmation  of  the  theory.  An  investigation  of  the  active 
loading  approach  described  here  when  used  against  realistic  scattering  bodies  (and  not  the 
surrogate  dipoles  described  here)  needs  to  be  conducted.  Such  an  investigation  could  be 
conducted  in  an  experimental  facility  or  with  numerical  electromagnetic  codes  with  the  objective 
of  comparing  the  measurements  made  with  the  probe  array  with  those  made  using  an  isolated 
probe  as  is  conventionally  done  in  near-field  bistatic  scattering  ranges.  The  investigation  would 
provide  confidence  in  the  accuracy  of  the  probe  array  approach  and  provide  data  to  estimate 
processing  requirements  and  time  reductions  implicit  in  the  probe  array  approach  for  near-field 
bistatic  measurements. 
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